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ABSTRACT 

Currently, the analysis of transmission spectra is the most successful technique to probe the 
chemical composition of exoplanet atmospheres. But the accuracy of these measurements is con¬ 
strained by observational limitations and the diversity of possible atmospheric compositions. Here 
we show the UV-VIS-IR transmission spectrum of Jupiter, as if it were a transiting exoplanet, 
obtained by observing one of its satellites, Ganymede, while passing through Jupiters shadow i.e., 
during a solar eclipse from Ganymede. The spectrum shows strong extinction due to the presence 
of clouds (aerosols) and haze in the atmosphere, and strong absorption features from CH 4 . More 
interestingly, the comparison with radiative transfer models reveals a spectral signature, which 
we attribute here to a Jupiter stratospheric layer of crystalline H 2 O ice. The atomic transitions 
of Na are also present. These results are relevant for the modeling and interpretation of giant 
transiting exoplanets. They also open a new technique to explore the atmospheric composition 
of the upper layers of Jupiter’s atmosphere. 

Subject headings: Solar System: Jupiter atmosphere - Exoplanets: characterization 


1. Introduction 


Over the past two decades, more than 1800 


exoplanets have been discovered, Borucki et al. 


(2010), approximately 65% of which are transit- 


^ Departamento de Astrofisica, Universidad de La La¬ 
guna, Av., Astroffsico Francisco Sanchez, s/n, E-38206 La 
Laguna, Spain 

^Fundacio Observatori Esteve Duran. Avda. Montseny 
46, Seva 08553, Spain 

^Institut de Ciencies de lEspai (CSIC-IEEC), Campus 
UAB, Facultat de Ciencies, Torre C5, parell, 2a pi., E- 
08193 Bellaterra, Spain 


ing. For a small sample of these planets - those 
which orbit bright stars, and have large planet- 
to-star area ratios - their atmospheres can be ex¬ 
plored through transmission spectroscopy. During 
the transit of a planet in front of a star the stellar 
flux is partially blocked, but a very small fraction 
of the stellar flux, 10“^ for a Jupiter-like planet 
and a sun-like star system, passes through the thin 
planetary atmosphere (if an atmospheric thickness 
of 0.1 times the radius of Jupiter is considered). 
Transmission spectroscopy has allowed the detec¬ 
tion of atmospheric Na I, H I, G II, O I, H 2 O, GH 4 , 
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and CO 2 ( 

Huitson et al.||2013 

Tinetti et al.||2010 

Bellucci et al. |20041 Borucki et al.||2010 

Ehrenre-| 

ich et al. 20061 Sing et al. 20111. These observa- 


tions, however, push the detection capabilities of 
space and ground-based observatories to the limit 
and the results are often a source of discrepan- 


cies in the literature ( 

Ehrenreich et al.||2006 Sing 

et al. 

2011 

Gibson et al. 

2011 

1 , which need to 


be addressed. Thus, observing planetary transits 
in our own Solar System can serve as an invalu¬ 
able benchmark, and provide crucial information 
for future exoplanet characterizations. 

Here, we report the transmission spectrum of 
Jupiter, with high signal-to-noise ratio, as if it 
were a transiting planet. Our technique is to 
observe Ganymede, which is in synchronous ro¬ 
tation around Jupiter, when crossing Jupiter’s 
shadow. During the eclipse, the spectral features 
of the Jovian atmosphere are imprinted in the sun¬ 
light that, after passing through Jupiters plan¬ 
etary limb, is reflected from Ganymede towards 
the Earth, see Figure The ratio spectrum of 
Ganymede before and during the eclipse removes 
the spectral features of the Sun, the local telluric 
atmosphere on top of the telescopes, and the spec¬ 
tral albedo of Ganymede. Ganymede and Europa 
are practically atmosphere-less bodies and do not 
introduce any significant variability in the spec¬ 
tra. Similar observations have previously been 
applied to retrieve the Earth’s transmission spec¬ 


trum, through lunar eclipse observations (Palle et 

ar]|2^ . 


2. Observations 


Initially, we observed an eclipse of Ganymede 
on 06/10/2012, using LIRIS (Manchado et al. 


20041 at WHT in La Palma Observatory, Spain. 
This experiment was later repeated, and the re¬ 
sults confirmed, by observing a second eclipse 
with XSHOOTER ( [Vernet et alT][2MT| ) at VLT in 
Paranal Observatory, on 18/11/2012. Here, we fo¬ 
cus the discussion on the VLT data, due to their 
higher signal to noise ratio, but a detailed analy¬ 
sis of the WHT observations leads us to virtually 
identical results. 


The larger aperture of VLT allowed us to take 
rapid measurements - these eclipses are only ob¬ 
servable from the ground typically a few of times 
per year and the suitable observing window last 


several minutes - with high spectral resolution. 
At the same time, we extended our measurements 
into the visible and the near-ultraviolet regions, 
covering from 300 to 2500 nm in a single expo¬ 
sure. This is possible through a dichroic split¬ 
ting of the beam in three arms: UV, VIS and 
near-IR. The 0.5”, 0.4” and 0.4” slit widths were 
used respectively for each range, providing av¬ 
eraged R = 9100, 17400 and 11300, respectively. 
The observing method consisted on taking unin¬ 
terruptedly spectra of Ganymede in stare mode 
during its translation through Jupiter’s shadow. 
The telescope active optics was used for guiding 
at Ganymede’s non-sidereal rate in order to keep 
the target properly centered within the slit width 
during the dark phases of the event. The apparent 
magnitude of Ganymede is V=4.6 but it decreases 
by about 8 magnitudes during the umbra phase 


Smith 

(19751 

Tsumura et al. 

(2014 


Observations started at 3:30 UT, with a see¬ 
ing below which remained stable during the 
night. The telescope flexure compensation pro¬ 
cedure was carried out at 04:00 UT in order to 
keep the three slits staring at the satellite during 
the darker phases. The penumbra phase started at 
04:19 UT (see Figure[^, in which the planet blocks 
part, but not all, the direct sunlight. During this 
phase, the amount of direct sunlight reaching the 
satellite decreases with the progress of the eclipse, 
whereas the light refracted by the planetary atmo¬ 
sphere reaching the satellite increases. The eclipse 
was total when the umbra phase started and all 
direct sunlight was blocked, at 04:40 UT. During 
the umbra, only refracted light from the plane¬ 
tary atmosphere reaches the satellite. The final 
out-of-the-eclipse phase was not observed due to 
the occultation of the satellite by Jupiter’s disk, 
which happened at 06:00 UTC. 

The data reduction of the XSHOOTER was 
performed using the ESO pipeline Reflex. Each 
individual exposure had enough signal to noise ra¬ 
tio to be analysed independently. Nevertheless, 
the three CCDs have different readout times, and 
we took the VIS exposures as references. We then 
averaged the several Near-IR spectra taken dur¬ 
ing a VIS exposure, and the closest UV spectra in 
time to construct a full spectrum from 300-2500 
p,m. This way, 41 individual spectra were deter¬ 
mined corresponding to sunlit Ganymede before 
the eclipse (Bl), 8 during the penumbra (PI), and 
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three during the Umbra (U), at all wavelengths 
simultaneously. 


The penumbra transmission spectra were deter¬ 
mined by calculating the ratio between the spectra 
taken during PI and the spectra taken just before 
the eclipse Bl. The umbra transmission spectra 
were determined by calculating the ratio between 
the spectra taken during U and the spectra taken 
just before the eclipse Bl. Again, these ratios will 
cancel the telluric contribution of the local atmo¬ 
sphere, the solar spectrum intrinsic in the observa¬ 
tions, and the spectral signatures of the satellite 


(Montanes-Rodriguez et al. (2006); Palle et al. 


(2009); Yang et al. (2014)), leaving only the con¬ 
tribution from Jupiter’s atmosphere. 


During these observations, the airmass de¬ 
creases as the eclipse progresses (from 1.565 in Bl 
to 1.510 in the deepest U), making it easier to iden¬ 
tify residual telluric features, because they would 
be seen in emission in the (pen-)umbra/bright 
ratios (see Figure [^. 


3. Model Simulations 


Transmission spectra of the occultation of the 
Sun through Jupiter’s atmosphere as observed 
from Ganymede have been computed, simulating 
the observations of the eclipse in the penumbra 
and within the first stages of the umbra. The 
extension (horizontal) of the Sun’s disk as it is 
setting on Jupiter’s horizon has been considered. 
Because of the strong refraction of Jupiter’s atmo¬ 
sphere, this makes it possible to sound Jupiter’s 
limb with moderate vertical resolution (a few tens 
of km) at the lowest tangent heights. 

The transmission spectra have been calculated 
by using the Karlsruhe Optimized and Precise 


Radiative Transfer Algorithm (KOPRA; Stiller, 


(2000)). KOPRA is a well-tested line-by-line ra¬ 
diative transfer model which offers all the neces¬ 
sary physics for studying this problem. This code 
was originally developed for use in Earth’s atmo¬ 
sphere and has been recently adapted to the atmo¬ 


spheres of Titan and Mars (Garcia-Comas et al. 


2011 Robert et al. 2012). Here, the reference 


Jupiter’s atmosphere, including pressure, temper¬ 
ature and the species abundances were taken from 
Gonzalez et ah, (2011). We include the major 
gaseous species, CH 4 and H 2 , and other minor as 
C 2 H 2 , CO, H 2 O and NH 3 . The concentration of 


NH 3 (only relevant in the lowest troposphere) was 
taken from [Griffith et a~ (1991 ).The only species 
with significant concentration that is not included 
is ethane (C 2 H 6 ). The absorption of its strongest 
band in the region of 2.2 — 2.5 /im would be, how¬ 
ever, masked by the stronger CH 4 bands in that 
region. The molecular spectroscopic data for all 
species have been taken from the HITRAN com¬ 
pilation, 2012 edition ( Rothman et ah, 112013 1. 

Rayleigh scattering by molecular hydrogen and 
helium have also been taken into account by in¬ 
cluding the Rayleigh optical cross sections pro¬ 


vided by Ford et al. (1973) for H 2 and by Chan 


et al. (1965) for Helium. Ro-vibrational absorp¬ 


tion bands resulting from collisions between pairs 
of H 2 -H 2 and H 2 -He, the so-called collisions in¬ 
duced absorption or CIA, are also significant in 
the lower Jupiter atmosphere, where they form a 
smooth feature mainly in the 2.0-2.5 /im region. 
Our simulations include this absorption with ab¬ 
sorption coefficients at low temperatures, derived 


by 

Borysow 

(2002 

) for H 2 -H 2 pairs, and bj/ 

Bo- 

rysow 

(1989 

and 

Borysow & Frommhold 

(1 

989) 


for H 2 -He. 

In addition, we also included in the simulations 
Mie scattering by water ice and aerosols. The op¬ 
tical properties of the crystalline water ice (real 
and imaginary part of the refractive index) were 
taken from |Mastrapa et al. (2008 1 For a tempera¬ 
ture of 150 K. The aerosol particles were assumed 
to have a mean radius of 0.25 /rm, in the range 
of 0.2-0.5 /rm derived by Zhang et al. (2013) for 
the equatorial particles. The real part of the re¬ 
fractive index of these aerosol particles was taken 


from 

Khare et al. 

(1984 

) and the imaginary part 

from 

Zhang et al. 

2013 



4. Results 

In Figurej^, the penumbra spectrum of Jupiter 
is shown, representing a direct transmission com¬ 
ponent through the upper layers of the atmo¬ 
sphere of the planet. This is directly compara¬ 
ble to the information retrieved through a plan¬ 
etary transit. As expected, the spectra show 
the major absorption features of the most abun¬ 
dant component, CH 4 . But most interesting is 
the change in the spectral continuum: at shorter 
wavelengths (UV-VIS) transmission is lower due 
to the extinction caused by aerosols. Extinction 
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by clouds and aerosol particles have a smooth 
wavelength-dependent extinction effect on the 
stellar flux which does not produce sharp features. 
Clouds and hazes have been tentatively detected 
in Hot Jupiter atmospheres using transmission 
spectroscopy from the Hubble Space Telescope 
Ehrenreich et al. (2006), Pont et al. (2012). Our 


results give additional confidence to these findings. 

In Figure [^, several umbra spectra of Jupiter 
are shown. As the eclipse progresses, the geome¬ 
try allow us to sample lower tangent heights, and 
thus probe lower into Jupiters atmosphere. The 
planetary atmosphere has also the effect of refract¬ 


ing the stellar light (Garcia-Muhoz et al. 2012). 


This refracted component becomes more impor¬ 
tant as the eclipse progresses deeper into the um¬ 
bra. Thus, these spectra present a much promi¬ 
nent absorption not only by the aerosol particles 
(Mie scattering) mainly in the near-IR, but also by 
the CH 4 absorbing bands. However, the most in¬ 
teresting feature is the possible detection of strato¬ 
spheric H 2 O ice cloud features at 1.5 and 2.0 pm. 

In Jupiter, the visible cloud deck is composed 
of NH 3 ice condensation, and spans from 700 
up to 200 mbar at the equatorial region, near 
the tropopause, which is especially conspicuous 
in planetary images taken in the 890 nm strong 
methane band ( West et al.||2004 ). Below the visi¬ 
ble cloud layer, clouds formed by NH 3 ice particles, 
Ni 74 -SH solid particles, H 2 O ice and H 2 ONH 3 liq¬ 
uid solution are also present, but they are hard to 


detect (Taylor et al. 2004). Above these levels, the 


major atmospheric constituent is CH 4 , but recent 
analysis of observations made with the Cassini 
Imaging Science Subsystem (ISS), indicates that 
the planet is wrapped up by a haze at the lower 
stratosphere at the equatorial region and mid lat¬ 
itudes, at « 50mbar, which rises at higher levels 


(p < 20m 6 ar) over the poles (Zhang et al. 2013). 
But the composition of this layer remained unde¬ 
termined. 

It is from these upper levels that our spectro¬ 
scopic signals are coming from. Over the umbra 
evolution, H 2 O ice absorption features at 1.5 and 
2.0 /im appear and later diminish (see Figure]^) 
as the sounded tangent height on the planetary 
limb crosses the altitude (or pressure levels) where 
these particles are located. Their detection in our 
transmission spectra allows us to determine that 
Jupiter’s upper atmospheric hazes near the equa¬ 


tor are at least partially composed of a cloud deck 
of H 2 O ice crystals. While these H 2 O ice absorp¬ 
tion bands are well known, the spectral measure¬ 
ments of VIMS on board of the Cassini spacecraft 


taken in 2000 did not detect them (Formisano et 

ar]|200^ . 


In Figure our model simulations of Jupiter’s 
umbra and penumbra transmission spectra are 
plotted. All major features of the measured 
spectra are simulated. The absorption spec¬ 
tral features at 1.5 and 2.0/rm can be very 
well reproduced in our model with crystalline 
water ice, needing a total column of about 
10^^particles/cm'^ with a size of 0.01 p,m located 
near the 0.5mbar level, a water amount 420 times 


larger than that measured by HERSCHEL, Cav- 


alie et al. (2013), which found a water mass load 
of 1.5xl0~'g/cm^ in the gas phase. While the 
good spectral match supports the identification 
of H 2 O ice, the large amounts and the location 
at low pressure implied by the data are puzzling. 
Moreover, the presence of water ice at 0.5 mbar 
does not have an obvious explanation. 

In terms of aerosol mass loading, |Zhang et al.| 


(|2013|) derived 1x10 ‘^gjcm'^ at low latitudes and 


lxl0~‘^g/cm^ at high latitudes. The H 2 O ice 
particle distribution we derive (log normal with 
sigma=0.3 and r=0.01/rm) gives a H 2 O ice mass 
load of 6.3x 10“®g/cTO^. Hence, although it is 
larger that their aerosol mass load for low lati¬ 
tudes, it is comparable to their mass load for high 
latitudes, using observations at 0.25 and 0.9 /rm. 
Hence, we do not have any conflict between the 
derived H 2 O ice load and the results from previ¬ 
ous literature. The haze parameters used in our 
model are in agreement to those recently reported 


for Jupiters equator by Zhang et al. (2013). 


The solar atomic lines should disappear in the 
ratio between the (pen- )umbra and bright spectra 
of Ganymede. However, some residual features re¬ 
main, as several effects will alter the shape of the 
spectral lines, including: i) small Doppler shifts 
due to the Sun-Ganymede-Earth relative speeds, 
ii) changes in the solar region contributing to the 
spectra, ii) instrumental flexures, and iv) Raman 
scattering producing a center-to-limb brightening 


Yang et al. (20141. This latter effect seems to dom¬ 


inate the spectral shape of the ratio spectra near 
atomic solar lines, where the line residuals have an 
inverted W-shape pattern. 
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In Figure]^ we plot several penumbra and um¬ 
bra spectra (both divided by the direct sunlit 
Ganymedes spectrum) focusing on the region near 
the Na I doublet. Most of the lines are strongly 
deformed in a classical W pattern, as is the case 
of the Fe lines in between the Na I doublet, and 
other solar lines at 586.2 and 591.4 /im. But in 
the case of the Na I, there is a net absorption in¬ 
dicating the presence of a Na layer in Jupiter’s 
upper atmosphere. This Na I absorption starts to 
be clearly detectable in the last penumbra spec¬ 
tra (at high SNR) and in all umbra spectra (with 
progressively decreasing SNR). We plan to carry 
out a detailed study of these features in a future 
paper. 


The presence of Na in Jupiter’s upper atmo¬ 
sphere can be explained by the deposition of ei- 


ther cometary impacts ( 

Noll et al. 

1995 

) or the 

continuous outward flux of Na from lo ( 

Hendillo 


sunlit Ganymede spectra do not show any signs of 
Na absorption, and neither do the first penumbra 
spectra, indicates that the origin of this absorp¬ 
tion is located within Jupiter’s atmosphere, and 
not in the torus of Na trailing along lo’s orbit. 


The WHT observations have not been discussed 
at lenght in this manuscript, however the results of 
that earlier campaign are similar to those obtained 
with VLT. In Figurethe transimission spectrum 
of Jupiter in one of the umbras of the WHT data 
is compared to one of the umbra obtained with 
the VLT. While there are differences in the spec¬ 
tra due to the the fact that the eclipse geometry 
is not exactly the same, and that WHT needs two 
different exposures (different times, different um¬ 
bra depths) to cover the 0.9 to 2.5 /rm range, they 
show essentially the same spectral signatures. 


atomic transitions of Na are also present. These 
results are relevant for the modeling and interpre¬ 
tation of giant transiting exoplanets, but they also 
open a new technique to characterize the upper 
layers of Jupiter’s atmosphere. Taking advantage 
of the scanning of faint features that this technique 
provides, observations of other satellite eclipses 
from space could help to set limits to the strato¬ 
spheric water abundance in the upper layers of 
Jupiter’s atmosphere and provide a way to moni¬ 


tor the rate of cometary impacts on Jupiter (Smith 


et al. 1999) which, in turn, has consequences for 


the formation history of the solar system. 
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Fig. 1.— Top: A not-to-scale diagram showing the orbital geometry of the Jovian System during the 
observations. These events are observable from the ground when the Earth is near quadrature. The ratio 
spectra of Ganymede when in the umbra (or penumbra) and when fully sunlit, results in the determination 
of the transmission spectrum of Jupiter. Bottom: Universal time (UTC) and Angular Separation from 
Ganymede center to Jupiter center (JAS) for the eclipse of the satellite on 18/11/2012, as seen from Paranal, 
Chile. 
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Fig. 2.— The transmission spectrum of Jupiter during the penumbra (panel a) and the umbra (panel b) 
phases. Gray shaded regions mark the deeply absorbed telluric bands of H 2 O, which cannot be observed 
from the ground. The locations of the major CH^ absorption bands are marked with a brown background 
shadow. On both panels, the thin gray line is the brightness ratio of the light reflected off Ganymede when 
fully illuminated by the Sun at two different airmass (equivalent to the airmass difference for the two spectra 
used to extract the transmission spectrum). It serves to illustrate the contamination from telluric lines to be 
expected during our measurements. All spectra have been binned to a lower resolution for display purposes. 
Panel a): Two of our penumbra spectra, taken at different airmass, which are almost identical except for 
the stronger effect of Jupiter’s haze absorption as the eclipse progresses. Panel b): Three measured umbra 
spectra, normalized to have unity flux at 1.6 /rm. The different relative depth of the bands is due to different 
timing within the eclipse umbra, which allows probing different altitudes in Jupiter’s upper atmosphere. The 
first two spectra (green) show features of gaseous C'iJ 4 and H 2 O ices, while in the third spectrum (red), the 
H 2 O ice features diminish indicating that we are probing below the H 2 O ice cloud layer. In the third, deep 
umbra spectra the spectrum is dominated by the CH 4 ^ absorption and by the haze. 
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Fig. 3.— Top) Transmission spectra of Jupiter calculated for the early phase of the penumbra, and at several 
stages during the umbra over the 0.5 — 2.5 /im spectral region. The spectra show the most prominent CH^ 
bands, the extinction of the aerosol particles (haze), and the two distinct absorption peaks of water ice at 1.5 
and 2.0 /rm. Plotted with a grey line is one of our observed penumbral spectra from Figure]^. Bottom) An 
illustration of the distribution of haze, CH^, and H 2 O ice in our model, together with geometrical information 
on the range of tangent altitudes covered in each of our simulations in the top panel, corresponding to the 
mean times when the (umbra)-penumbral spectra were taken. Altitude is taken as zero at the 1 bar pressure 
level. In addition, we have over-plotted the profiles of the concentrations of aerosols and of the water ice 
particles cloud used in the computation of the transmission spectra. Note how the sounded region shrinks 
as the eclipse progresses, while moving down to lower tangent heights. 
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Fig. 4.— The evolution of Ganymede spectra during the observations centered on the Na I doublet wavelength 
region. The grey spectra are fully illuminated Ganymede spectra immediately before the penumbra phase, 
divided by a previous spectra also in full sunlight illumination (arbitary reference spectrum). These spectra 
are basically featureless. In blue are the penumbra spectra divided by the same reference spectrum, and 
they start to show signs of extra absorption at 0.5889 and 0.5896 /rm. In green are the two umbra spectra 
for which the SNR is good in this wavelength region. 
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Fig. 5.— The transmission spectrum of Jupiter as measured with the WHT telescope (black line) during 
the umbral phase. Gray shaded regions mark the deeply absorbed telluric bands of H 2 O, which cannot be 
observed from the ground. As in Figure the thin gray line is the brightness ratio of the light reflected 
off Ganymede when fully illuminated by the Sun at two different airmass. Overplotted (red line) is one of 
Jupiter’s umbral spectra taken with the VLT. Note how both WHT and VLT spectra shown the fingerprints 
of gaseous CH^ and H 2 O ices. 
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